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Study on the specific heat of wood by statistical mechanics 
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Abstract: From the microstructure of wood, theoretical expressions of the wood specific heat were derived by 
statistical mechanics. With the theoretical expressions derived, the theoretical values of specific heat for 33 tree 
species, with different moisture contents and under varied temperature conditions were calculated and compari¬ 
son was also made to the experimental values under the same conditions. The results showed that the maximum 
error and mean error by the theoretical expressions of this paper are only 7.8% and 2.5% respectively, while 
those error of the theoretical values for 33 tree species calculated with Dunlap's empiric equation were 15.2% 
(max.) and 9.3% (mean), and forKHpH^^OB empiric equation, they were 20% (max.) and 11% (mean). 
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Introduction 

Since wood is a natural macromolecular organism, 
and chemical composition, structure and biological 
cell structure of wood are extremely complicated, 
these cause great difficulty to the theoretical study on 
the specific heat of wood. So far, any research 
achievement has not been reported on it. 

We try theoretically to find out the internal relation 
and essential laws between the specific heat of wood 
and the microstructure of wood, reveal the internal 
factors and the external conditions of affecting the 
specific heat of wood, and finally derive a theoretical 
expression of specific heat of wood, it can be used to 
replace the experimental which measurement of the 
wood specific heat. 


volume of each glucose radical is 

v = — ahcstn 3 =0.1689X 10’ 27 m 3 
4 

Theoretical expression of specific heat of 
wood 

First let’s consider such a molecular crystal which 
contains N molecules and each molecule is com¬ 
posed of S atoms. It is assumed that there are 3 N 
modes of vibration whose frequency distribution can 
be approximated according to Debye’s theory and 
the frequency distribution of the other 3A/(S-1) modes 
of vibration can be approximated according to Ein¬ 
stein’s theory. Then the logarithm of the crystal parti¬ 
tion function K(T) can be expressed as 


Microstructure of wood 

The physical property of wood depends mainly on 
substances composing the cell wall of wood, and 
cellulose is the major composition of the cell wall. 
Proved by X-ray diffraction, all natural cellulose have 
the same X-ray diffraction pattern, which indicates 
that all natural celluloses have the same crystal 
structure. The smallest unit of the crystal structure of 
cellulose is cell. The model of unit cell belongs to the 
monoclinic system, and the length of axis is a=0.835 
nm, b=1.03 nm, c=0.79 nm respectively and the an¬ 
gle p between a axis and c axis is 84°. There are four 
glucose radicals (CgHioOs) in each cell, and the 
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Where E 0 is zero point vibrational energy of the 
crystal. From expression (1) we can get the heat ca¬ 
pacity of the molecular crystal at constant volume. 
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Where 0 = hv > _q ) _ p( & r> ) are respec- 

k T 

tively called Debye characteristic temperature and 
Debye function. Their expressions are respectively 
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(4) 


Where u is sound velocity of wood, v is each 

molecules’volume and x = — . 

kT 

But the problem is that the celluloses can’t be sim¬ 
ply regarded as the crystal composed of the mole¬ 
cules, for every glucose radical doesn’t lie in the 
similar place of the lattice. On the average, every 
four glucose radicals forms a cell which repeats itself, 
namely the cellulose belongs to the complicated lat¬ 
tice crystal. So we can’t make direct use of expres¬ 
sion (2) to calculate the specific heat of wood. 

In Born’s view (Qian 1962), a crystal containing N 
cells with SN atoms can have 3s groups of different 
normal modes of vibrations and frequencies. The 
wave numbers of each group is uniform distribution 
in the three-dimensional space. To every wave 
number, there is correspond to two frequencies, a 
higher frequency which hardly depends on the wave 
number to remain a constant is called optical mode 


of vibration, a lower frequency which almost keeps 
direct proportional to wave number is called acousti¬ 
cal mode of vibration. In the crystal made up of N 
single cells, each group composed of N acoustical 
modes of vibration contributes a Debye term, while 
each group composed of N optical modes of vibra¬ 
tion contributes an Einstein term. So any crystal par¬ 
tition function can roughly be expressed with the 
product of 3S partition functions. Suppose in 3S par¬ 
tition functions there is 3a belong to the Debye type 
(Here a is equal to one at least) and all the other 
3(s-a) belongs to the Einstein type. For this rea¬ 
son expression (2) can be changed into 


C v =3 aNkD(^-) + Nk 


3 s 


i=3a +1 


(Qj/IT) 2 
sinh 2 (0,- I2T) 


(5) 


The determination of how many groups of vibration 
modes in the 3S groups belong to the acoustical 
mode and how many groups belong to the optical 
mode depends on the particular structure of the 
crystal. As far as the cellulose crystal of wood is 
concerned, it is linear macromolecule. The basic unit 
of the cellulose structure is D-glucose. The cellulose 
is a chain-shaped molecule joined by P -1, 
4-glucosidic linkages without any branches. The 
chemical structure of the big cellulose molecule is 
represented in Fig. 1. 


ch 2 oh 




ch 2 oh 






Fig. 1 Chemical structure of the cellulose 


The cellulose molecule in the long chain direction 
combines to form a succession structure depending 
on valence bonds, and in side direction it forms as¬ 
sociated molecular bundle depending on hydrogen 
bonds. It is the parallel arrangement of the molecular 
chain with excellent direction and the associated hy¬ 
drogen bonds in the side direction that form crystal 
structure of wood cellulose. 

According to the above wood cellulose’s chemical 
constitution and crystal structure, we can put forward 
the microstructure model of cellulose crystal to derive 
specific heat of wood, as shown in Fig.2. It is clear 
that a pyran-ring is made up of 20 atoms, 60 degrees 
of freedom, 54 degree of freedom of internal vibra¬ 
tion, including the stretching and bending vibration of 


C—C, C—H, C—O, O—H bonds. They all belong to 
the optical mode of vibration, and each contributes 
an Einstein term. As to the other 6 freedom degrees, 
one half are translation degree of freedom in the 
mass center of pyran-ring, the other are rotary de¬ 
gree of freedom. Because the pyran-ring is a looped 
structure as and too large to be considered as a par¬ 
ticle, we have to consider its rotary motion. These 6 
degrees of freedom satisfy the dispersion relation of 
acoustical motion and each contributes a Debye term. 
As to the oxygen atoms between the pyran-rings, the 
intensity of C—O bond is similar to that inside the 
pyran-ring, but the bonds on the pyran-ring form a 
cyclic construction which can be considered as a 
whole unit. At the same time, the mass of the oxygen 
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atom between the two pyran-rings is much smaller 
than the total mass of the atoms in the pyran-ring. 
Therefore, it can’t form Debye term. 



O: oxygen atom-hydrogen bond AlW glucosidic bond 


where e 


f 0 O 
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/ sinh 


{ T J 


2 T/ 


[it jj 


is Einstein’s func¬ 


tion. Substituting known quantity into Eq.(8), we still 
further obtain more simply theoretical expression of 
the wood specific heat 


C 0 (/) = l.2294 +0.006714/ (kJ- kg' 1 - 1C 1 ) (9) 

The relationship between the specific heat of the 
wet wood C w (whose moisture content is w%) and 
that of the absolutely dry wood C 0 (it is expressed in 
Equation (9)) is given by the equation 


„ C 0 + C'W% 

w ~ i + w% 


( 10 ) 


Fig. 2 The model of cellulose structure for Comparison between the theoretical value 

theoretically deriving wood specific heat experimental value 


So these 3 freedom degrees of the oxygen atom 
between the two pyran-rings belong to the optical 
mode of vibration and each contributes an Einstein 
term too. Thus, equation (5) can be turn into 


C v = 




v=\ 


(Q V /2T) 2 

sinh 2 (0 u /27 r ) 


( 6 ) 


We calculated the specific heat of 33 species of 
trees with different moisture contents and tempera¬ 
ture by means of equations (9) and (10), and made a 
comparison between the theoretical values of spe¬ 
cific heat and their experimental values under the 
same condition (temperature and moisture content). 
The results are shown in Table 1. 

Conclusion 


In this way we can obtain the theoretical expres¬ 
sion of dry wood specific heat. 



{Qy/2 T) 2 
sin b 2 (Q v /2T) 


(7) 


Where R is universal constant of gases, M is 
molar mass of a glucose radical, and 0 U is the 

characteristic temperature of each chemical bond in 
a pyran-ring. 

In order to simplify calculation, we can expand ex¬ 
pression (7) into Taylor series at T=T 0 (T 0 =273.15 K). 
When the series is kept to the first order of solution, 
we obtain 


The comparison between the theoretical values 
and the experimental ones of the specific heat in the 
33 species of wood with different moisture contents 
and temperature in Table 1 shows that the maximum 
error is 7.8%, and the average error is 2.5%. 

If the specific heat of the 33 species of wood in Table 
1 is calculated by Dunlap’s (Cheng et at. 1985) em¬ 
piric equation 

Co = 0.266 + 0.00116/ (cal- g" 1 K' n ) 


Or by Kh pH/in ob empiric equation 


C 0 = 0.28 

v 



0.2 


+ 0.09 


(cal- g- 1 K- ') 



it can be found that the maximum error of the former 
is 15.2% and that of the fatter 20%, and the average 
error of the former is 9.3% and that of the latter 11%. 

As shown above, theoretical expression of wood 
specific heat that is derived using statistical thermo¬ 
dynamics from the microstructure of wood is a more 
accurate formula. 
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Table 1. Theoretical and experimental values of specific heat in 33 species of wood 


Species of trees 

Moisture content 

/% 

Temperature 

Specific heat /kJ 

■ kg' 1 • K' 

Error 

/% 

rc 

Theoretical 

Experimental. 

Picea jezoensis var. komarovii 

12.7 

16.5 

1.661 

1.624 

+2.3 

Pinus koraiensis 

11.0 

16.5 

1.623 

1.666 

-2.6 

Cunninghamia lanceolata 

17.0 

23.5 

1.794 

1.762 

+1.8 

Pinus massoniana 

14.6 

24.2 

1.748 

1.737 

+0.5 

Podtocarpus imbricatus 

13.6 

17.4 

1.687 

1.687 

0 

Larix olgensis 

14.5 

13.9 

1.686 

1.720 

-2.0 

Paulownia fargesii 

10.5 

17.2 

1.615 

1.670 

-3.3 

P fortunei 

13.0 

21.7 

1.699 

1.725 

-1.5 

P. elongatei 

12.4 

22.6 

1.691 

1.754 

-3.6 

P. austrais 

12.9 

21.2 

1.693 

1.649 

+2.7 

P kawakamii 

11.0 

18.7 

1.636 

1.649 

-0.8 

P. catalpifolia 

12.7 

22.4 

1.696 

1.770 

-4.2 

P. tomentosa var. tsinlingensis 

11.3 

17.6 

1.636 

1.695 

-3.5 

P tomentosa 

13.0 

22.4 

1.703 

1.712 

-0.5 

Populus euramericana 

12.1 

18.3 

1.658 

1.633 

+ 1.5 

Gossampinus malabaricus 

15.1 

23.7 

1.755 

1.628 

+7.8 

Alniphyllum fortunei 

13.9 

24.5 

1.734 

1.725 

+0.5 

Tilia mandshurica 

10.9 

16.0 

1.617 

1.615 

+0.1 

Juglans mandshurica 

14.5 

16.1 

1.699 

1.720 

-1.2 

Catalpa duolouxii 

17.3 

26.0 

1.814 

1.875 

-3.2 

Populus simonii 

13.3 

15.8 

1.671 

1.636 

+2.1 

Melia azedarach 

15.3 

24.5 

1.764 

1.720 

+2.6 

Manglietia hainanenss 

12.7 

21.0 

1.688 

1.712 

-1.4 

Ulmus sp. 

10.1 

25.0 

1.653 

1.670 

- 1.0 

Betula platyphylla 

13.2 

16.7 

1.674 

1.633 

+ 2.5 

Acer davidii 

12.3 

24.5 

1.699 

1.716 

-1.0 

Acer mono 

11.7 

25.0 

1.689 

1.712 

-1.3 

B. alnoides 

15.3 

24.2 

1.762 

1.746 

+0.9 

Fraxinus mandshurica 

12.1 

26.0 

1.704 

1.771 

-3.8 

Quereus mongolica 

12.7 

25.0 

1.711 

1.846 

-7.3 

Schima wallichii 

16.0 

23.6 

1.774 

1.725 

+2.8 

Homalium hainanense 

12.0 

24.0 

1.690 

1.771 

-4.6 

Quereus acutissima 

14.8 

21.0 

1.733 

1.880 

-7.8 
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